Atomtronic circuits of diodes and transistors 
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We illustrate that open quantum systems composed of neutral, ultracold atoms in one-dimensional 
optical lattices can exhibit behavior analogous to semiconductor electronic circuits. A correspon- 
dence is demonstrated for bosonic atoms, and the experimental requirements to realize these devices 
are established. The analysis follows from a derivation of a quantum master equation for this general 
class of open quantum systems. 
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Atomtronics [1] focuses on ultracold atom analogs of 
electronic circuits and devices. Previous work has es- 
tablished the possibility of diode 0, Q and transistor 
0)0,0] behavior in atomic systems. We calculate the de- 
vice characteristics of these semiconductor-like systems 
by construction of complete atomtronic circuits. Such 
circuits contain the analog of electronic power supplies 
or batteries, and the necessary device connections. This 
opens up the possibility for more complex circuits in 
which atomtronic device components can be cascaded. 
The approach we introduce develops a general theoretical 
method for solving the dynamics of strongly-interacting, 
many-body, open quantum systems. 

A diode is in essence a device which exhibits an asym- 
metric response; applying a potential gradient one way 
leads to a large current and the opposite way leads to 
virtually no current. The atomtronic diode we present 
is novel compared to some previously proposed devices 
[3, 01 in that it does this without relying on an intrinsic 
irreversibility of the device itself {e.g. the spontaneous 
emission of a photon or phonon). The atomtronic transis- 
tor presented here contrasts with other proposals [1, Q , 
in that it operates in steady-state, rather than only in an 
initial (transient) regime. It also reproduces the essen- 
tial transistor behavior, including the control of a larger 
atomtronic current with a smaller one, the ability to per- 
form digital logic, and exhibits linear amplification. One 
motivation for designing atomtronic diode and transistor 
components that are intrinsically reversible is that they 
can then be combined in quantum computing applica- 
tions where coherent logic is required. 

This topic is relevant for the emerging field of strongly- 
correlated dynamical and nonequilibrium phenomena in 
ultracold atomic gases 0, 0, 0, Atomtronic de- 

vices require site- by-site control of the temporal or spatial 
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a common 



properties of optical lattices [111 
goal for numerous experiments aiming to utilize optical 
lattices for quantum control, quantum transport, quan- 



tum computing, and quantum simulation 15|, Il6l . Il7 |. 
What is novel here is that the atomtronic systems we 
consider are open quantum systems, meaning that they 
are composed of optical lattices coupled to two or three 



reservoirs, with the reservoirs acting as sources or sinks 
for particles. 

The Hamiltonian takes the separable form 



H = H„ 



Hrcs + Hv, 



(1) 



where Hgys and 



denote system and reservoir Hamil- 
tonians, respectively, and Hy their interaction. Although 
atomtronic devices can utilize either fermionic or bosonic 
atoms, we focus our discussion on the bosonic case, pro- 
viding contrast with electronics where the carriers are 
fermionic. In the lowest band approximation, TJgys is 
the Bose-Hubbard Hamiltonian that describes ultracold 



bosons in a one-dimensional optical lattice |12l . Il8l . \li 

given by 



ii,3) i 



u 



where site i has energy e^, interaction energy U , annihi- 
lation and number operators di and Ni = a\di respec- 
tively, and Jij is the hopping energy between adjacent 
sites i and j. To achieve a steady-state current, atom- 
tronic circuits require two or more reservoirs held at dif- 
ferent chemical potentials. The free Hamiltonian for the 
reservoirs is 



(3) 



where v identifies the mode of reservoir I with energy huj^i 
and annihilation operator R^i. It is assumed that each 
reservoir is so large that its thermodynamic properties 
are parametrized by a constant chemical potential /i; and 
temperature Tj. Each reservoir is connected to a single 
system site si so that the interaction between the system 
and reservoir can be written as 



Hv ^^guiRli^si +h.( 



(4) 



where g^i are the system-reservoir hopping energies. 

Our atomtronic analog of electronic circuits begins 
with the definition of the atomtronic battery, composed 




Reservoir {Mott insulator 
with unit filling) 



Reservoir (Mott insulator 
with zero filling) 



FIG. 1: (Color online.) (a) An electronic circuit with a bat- 
tery, resistive wires, and a p-n junction diode in forward bias, 
and (b) the atomtronic equivalent. The left reservoir has 
chemical potential fiL and unit filling, and supplies atoms 
that flow through the system (dashed box) to the right reser- 
voir, with chemical potential and zero filling. The weak 
system-reservoir coupling is ensured by large tunneling barri- 
ers at each interface. 



of two reservoirs of ultracold atoms having different 
chemical potentials, corresponding to different electric 
potentials at the terminals of a conventional battery. 
In Fig. [IJa), we show an electronic circuit composed of 
a diode connected to a battery by resistive wires. In 
Fig. [Hb) , the atomtronic equivalent of this diode is the 
part of the lattice that is enclosed by a dashed box with 
the lattice site energy on the left half offset from the lat- 
tice site energy on the right half. To the left and right of 
the diode are the reservoirs at high and low chemical po- 
tentials respectively. These are the analogs of the battery 
terminals. The system-reservoir separation is provided 
by large tunneling barriers that ensure weak coupling. 
The explicit physical implementation for the reservoirs 
shown is an optical lattice in the zero-temperature Mott- 
insulator phase. In this system, strong interactions lead 
to a fermionization of the bosons giving a broad contin- 
uum of reservoir modes fully occupied below the chemical 
potential. 

In the situation described we find that the current re- 
sponse depends strongly on the direction of the applied 
chemical potential difference. The atomtronic diode's 
conduction asymmetry becomes optimal when the height 
of the potential step in the middle of the diode is equal to 
the on-site interaction energy U. This general behavior 
is observed not only for the six-site lattice shown, but 
also for lattices containing variable numbers of sites. In 
order to explain the general mechanism, we consider the 
simplest case of a two- well system with potential offset U 
and connected on each side to a reservoir. An appropri- 
ate basis for this two-site case is the number basis |n, to), 
where n and to are the numbers of atoms in the left and 
right wells respectively. 



FIG. 2: (Color online.) Dynamical path of an initially empty 
two-well lattice system. The vertical axis corresponds to the 
energy of the system's quantum states. The solid (blue) and 
dashed (red) arrows indicate the transfer of atoms from the 
left and right reservoirs respectively, (a) Reverse bias (low 
left hand chemical potential and high right hand chemical 
potential). Here the system evolves into the |0, 2) dark state 
where it becomes trapped, (b) Forward bias (high left hand 
chemical potential and low right hand chemical potential). 
The resonance between |2, 0) and |1, 1), shown by the double- 
headed arrow (green), allows the resonant transport of atoms 
across the lattice. 



As shown pictorially in Fig. [2l the potential offset sets 
up a degeneracy between two states. The state |2,0) is 
resonantly coupled to the state |1, 1) (double headed ar- 
row), and both remain separated in energy from |0,2). 
The reverse bias situation is shown in Fig. [D^a). If the 
chemical potentials for both the left reservoir, fi^^ and 
right reservoir, fin, are initially sufficiently low, the sys- 
tem begins in an equilibrium devoid of particles. Raising 
/ifl from this point leads to populating the |0, 1) state. 
Since this state is detuned from |1, 0) by the on-site inter- 
action energy U, the steady state dynamics corresponds 
to virtually no current flow. Further raising fin allows 
the state |0, 2) to be occupied. This state is detuned 
from |1, 1) by 2U. As a result, there is again no effective 
transport across the lattice system. 

However, if fifj is reset to its initial low value and fiL 
raised, the resulting forward bias dynamics arc quite dif- 
ferent. As shown in Fig. (Hb), when fi^ is increased to 
the point where the |2,0) state is occupied, this opens 
up the possibility for resonant transport of atoms from 
the left reservoir to the right reservoir via the two res- 
onant current loops; |2,0) |1,1) |1,0) |2,0) and 
|2,0)^|1,1)^|2,1)-.|2,0). 

This qualitative picture of the diode operation is sup- 
ported in detail by the solution of the quantum dynamics 
for the open quantum system. We solve for the evo- 
lution of the reduced density operator of the system, 
a = Trrcs[p], in which p is the full density operator of 
the system and reservoir, and the trace is over the reser- 
voir degrees of freedom. Under the Born-Markov ap- 
proximation, the evolution of a is given by the master 
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equation (l^, [HI, [111 , 



+ ihsys j a{t) ^ J dT^{T)&{t), 



where the memory kernel is 

J'(t) = Trres i^V exp[-iT(Lsys + Lres)]«Ly 



(5) 



Here the Liouville operators Lgys, Lres, and Ly are de- 
fined via commutators of their respective Hamiltonians, 
i.e. for operator 6, KLxQ = [-ffjc, 6], X G {sys, res, V}. 



Substituting Eq. ([6]) into Eq. ([5]) yields the following 
closed form that is sufficiently simple to solve numeri- 



(6) cally: 



da 



1 

ih 



sys, < 



I 



(7) 



with as,(— r) = exp(— ii/sysT)asjexp(i_ffsysT) and 
Rvi{~t) = exp{-iHrcsT)R„icxp{iHT-cBT) defining the sys- 
tem and reservoir operators respectively in the interac- 
tion picture. The parameter rj regularizes the integral 
and remedies the infrared divergence that is needed to 
account correctly for the beyond Born-Markov effects 
that arise from the coupling of near degenerate system- 
reservoir levels [23|. The theory should be cutoff and 
renormalized to account properly for the ultraviolet di- 
vergence that arises from the replacement in Eq. ([7]) of 
g^i by one that is mode independent, i.e. g^i — gi. 

The numerical solution of the diode connected to a 
reservoir of this type in steady state is shown in Fig. [31 
The reverse bias (inset) shows virtually no current re- 
sponse as a function of the chemical potential difference, 
and the forward bias case exhibits a turn-on of the cur- 
rent at the point where the resonance states previously 
discussed become occupied. As anticipated from the 
qualitative picture given earlier, when the current turns 
on, there are two observable steps associated with linear 
combinations of the two possible resonant current loops, 
with a rise broadened by 77. 

Building on this understanding of the diode we now 
develop the atomtronic transistor. Since a transistor is a 
three terminal device, transistor action requires at least 
three lattice sites connected to three independent reser- 
voirs. As seen in Fig. IH the middle well is shifted down 
in energy by U with respect to both the left and right 
wells. To investigate the device characteristics, we put 
a fixed bias across the transistor. This means that /i^ 
is higher than the value necessary to give an occupancy 
of one atom on the left system site, and /i^ is below the 
value needed to deplete all atoms on the right system 
site. We then study the system's steady-state response 
to an increase of the middle chemical potential, ■ 

When /iAf is low, the middle site contains virtually no 
atoms and the system is locked in the |1, 0, 0). This state 
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FIG. 3: (Color online.) Current of the forward-bias and 
reverse-bias (inset) response for the diode. The horizontal 
axes are the chemical potential difference normalized to the 
resonance value. The vertical axes are the resulting current 
responses normalized to the saturation current 7s that arises 
when the normalized chemical potential difference is unity. 



is only weakly coupled to |0, 0, 1) via a second-order off- 
resonant process through |0, 1,0). As the middle chem- 
ical potential is raised, a non-zero occupancy results on 
the middle site, and a resonance is accessed among the 
states |1,1,0), |0,2,0) and |0, 1,1). This triggers a cy- 
cle that corresponds to a significant current across the 
system, as shown in Fig. [H High linear differential gain 
(inset) is achieved by ensuring that the system-reservoir 
coupling gi in Eq. ([7]) is smaller in magnitude for the 
middle reservoir compared to the left and right couplings. 
The current that results from the |0, 2, 0) to |0, 1, 0) tran- 
sition, i.e. an atom leaving the system and entering the 
middle reservoir, is then much smaller than the current 
which flows from left to right across the device. In Fig.[5l 
results are plotted of the numerical simulation for the 
case in which the middle reservoir coupling is one-fifth 
that of the reservoirs on either end. 

In conclusion, we have calculated device characteris- 
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FIG. 4: (Color online.) The level scheme for the atomtronic 
transistor with a fixed chemical potential difference across the 
system. The base potential is set to put one particle on the 
middle site, which triggers a resonance that allows for trans- 
port across the system. The left, right, and middle reservoir 
dynamics are represented by solid (blue), dashed (red), and 
dotted (purple), arrows respectively, while intrasystem dy- 
namics are represented by double-headed arrows (green). 
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FIG. 5: (Color online). The current response and differential 
gain of the atomtronic transistor are plotted. In this simula- 
tion, the reservoir couplings gi for the left and right reservoirs 
are five times that of the middle reservoir. For fixed Mz, — ^J^R, 
we vary fj,M and record the response of currents leaving the 
system from the right site (blue solid curve) and from the mid- 
dle site (green dashed curve). Here, 7s is the current value 
of the emitter when fj,M = ^m- The inset is a plot of right 
current versus the middle current, from which a large linear 
differential gain is apparent. 

tics of simple diode and transistor atomtronic circuits. 
We have focused here on the two-site and three-site mod- 
els to simplify the presentation, but full calculations for 
more extended systems with a greater site-number con- 
firm that the diode-like and transistor-like behavior per- 
sists. This work motivates expansion of the analogy to 
more complex circuits, such as amplifiers, constant cur- 
rent sources, flip-flops, and logic elements. Apart from 
their potential practical value, it would be intriguing to 
consider the connection of atomtronics with quantum in- 
formation physics and quantum computation, since com- 
puters are based on cascaded logic devices and these 
atom devices presented operate exclusively in the coher- 



ent regime. 

Atomtronic devices will always be much slower than 
their electronic counterparts. For example, the overall 
current flow is determined principally by the tunneling 
rate from one site to the next Compensating for 
this is the fact that atomtronic studies are pursuing de- 
vice physics in a novel physical system: unlike electrons, 
atoms have a complex internal structure and internal 
states, the atoms can be bosons or fermions, they are 
massive and affected by gravitational and electromag- 
netic fields, the material itself which is the optical lattice 
can be dynamically and spatially varied, and many un- 
avoidable effects present in solid-state systems such as 
crystalline impurities, dislocations, and phonon scatter- 
ing are absent. 
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